MPI-PhT/2000-45 
AZPH-TH/00-03 
November 2000 



> 
On 



X 



Percolation and the existence of a soft 
phase in the classical Heisenberg model 



o 
o 
o 

> 
o . 

^ ■ Adrian Patrascioiu 

(N 
(N 



Physics Department, University of Arizona 
Tucson, AZ 85721, U.S.A. 

and 

Erhard Seiler 

O . 

Max-Planck-Institut fur Physik 
- Werner-Heisenberg-Institut - 
Fohringer Ring 6, 80805 Munich, Germany 



Abstract 

We present the results of a numerical investigation of percolation prop- 
erties in a version of the classical Heisenberg model. In particular we 
study the percolation properties of the subsets of the lattice correspond- 
ing to equatorial strips of the target manifold S 2 . As shown by us several 
years ago, this is relevant for the existence of a massless phase of the 
model. Our investigation yields strong evidence that such a massless 
phase does indeed exits. It is further shown that this result implies 
lack of asymptotic freedom in the massive continuum limit. A heuristic 
estimate of the transition temperature is given which is consistent with 
the numerical data. 



PACS: 64.60. Cn, 05.50.+q, 75.10.Hk 



1 Introduction 



If one looks at textbooks to learn about the phase diagram of the two dimen- 
sional (2D) O(N) models the situation seems clear: for N = 2, there is a 
transition to a low temperature phase with only power law decay of correla- 
tions, whereas for the nonabelian case N > 2 there is exponential decay at all 
temperatures. But while the first statement has been proven rigorously a long 
time ago the second one remains an open mathematical question 0. The 
standard belief is rooted in the perturbative asymptotic freedom of the models 
for N > 2; but over the years we have brought forth many reasons why we 
think it is unfounded []3], f|, |5| . The absence of a mathematical proof together 
with ambiguous numerical results left the issue wide open. 

In this paper we would like to present what we regard as convincing nu- 
merical evidence that in fact the 2D 0(3) model possesses a massless phase 
for sufficiently large (3 and give a rigorous proof that this is incompatible with 
asymptotic freedom in the massive phase. We will also give a heuristic expla- 
nation of why and where the phase transition happens. 

The models we are considering consist of classical spins s taking values 
on the unit sphere S N ~ l , placed at the sites of a 2D regular lattice. These 
spins interact ferromagnetically with their nearest neighbors. Let (ij) denote 
a pair of neighboring sites. We will consider two types of interactions between 
neighbouring spins: 

• Standard action (s.a.): = —s(i) ■ s(j) 

• Constrained action (c.a.): = —s(i) ■ s(j) for s(i) ■ s(j) > c and 
Hij = oo for s(i) ■ s(j) < c for some c G [—1, 1). 

The corresponding Gibbs measures are (for a finite lattice) given by 

^.a. = |n^n^(^)) (i) 

(ij) « 

for the standard action and 

d»c. a . = \ n wo • so - ) - c)] n du ( a (0) (2) 

(ij) i 

for the constrained action, where di> is the standard measure on the two sphere 
S* 2 and the product Yluj) is over nearest neighbors. 

Almost a decade ago we showed || that one can rephrase the question of 
the existence of a soft phase in these models as a percolation problem and in 
fact this is the reason we introduced the c.a. model. It should be noted that 
the c.a. model shares with the s.a. model not only invariance under 0(N), but 
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has also the same perturbative (= low temperature) expansion and the same 
'smooth' classical solutions as the s.a. model. It is therefore to be expected 
that the s.a. and c.a. models fall in the same universality class (possess the 
same continuum limit) and, as we shall show shortly, the numerical evidence 
supports this expectation. The advantage of studying the c.a. model stems 



from the following fact: let e c = y 2(1 — c) and S ec the set of sites such that 
\s ■ n\ < e c /2 for some given unit vector n. Our rigorous result || was that 
if for some e > e c the set S e on the triangular (T) lattice does not contain a 
percolating cluster, then the O(N) model must be massless at that c. For the 
abelian 0(2) model we could prove the absence of percolation of this equatorial 
set S ec for c sufficiently large || (modulo certain technical assumptions which 
were later eliminated by M. Aizenman |7]]). For the nonabelian cases we could 
not give a rigorous proof. We did however present certain arguments ||, |9| 
explaining why the percolating scenario seemed unlikely. 

In this paper we will present numerical evidence that there is an e) such 
that for e < to S t does not percolate for any c; for sufficiently large c this eo 
will be larger than e c and the model will thus be massless. We will also show 
that due to a rigorous inequality derived by us in the past [[UJ , the existence of 
a finite (3 crt in the s.a. model on the square (S) lattice is incompatible with the 
presence of asymptotic freedom in the massive continuum limit of the model. 



2 Percolation and masslessness 

In this section we briefly review the special features of percolation in two 
dimensions and give a brief sketch of our argument relating percolation prop- 
erties to the absence of a mass gap. We restrict the discussion to the T lattice; 
this keeps the arguments simpler because the T lattice is self-matching and no 
distinction has to be made between connectedness and ^-connectedness (where 
points are also considered connected along diagonals). 

The following two facts special to 2D are relevant for our discussion: 

1. Noncoexistence of disjoint percolating sets: Let A be the subset of the 
lattice defined by the spin lying in some subset A C S N ~ l and A its comple- 
ment. Then with probability 1 A and A do not percolate at the same time. 
This has been proven rigorously only for special cases like Bernoulli percola- 
tion and the + and — clusters of the Ising model, but is believed to hold quite 
generally. (Aizenman J7[ showed that in the case of 0(2) one does not need to 
invoke this principle). 



2. Russo's lemma pl| : If neither A nor its complement A percolate, then 
the expected size of the cluster of A attached to the origin, denoted by (A), 
diverges; the same holds for its complement A. (In this simple form the lemma 
only holds for a self matching lattice like the T lattice). If A percolates, then 
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(A) is expected to be finite. 

The subsets of the sphere S 2 that interest us here are the following: 

• 'equatorial strip' <S e , defined by \s-n\ < e/2 for some fixed unit vector n. 

• 'upper polar cap' P+, defined by s ■ n > e/2, 

• 'lower polar cap' V~ , defined by s ■ n < —e/2. 

• 'union of polar caps' V e = U V~ . 

The subsets of the lattice defined by these subsets of the sphere we denote 
by the corresponding roman letters S e etc. and for brevity we say 'a certain 
subset of the sphere percolates' instead of 'the subset of the lattice induced by 
a certain subset of the sphere percolates' etc.. 

According to the discussion above, there are the following possibilities: 
either S t percolates, or P e percolates, or neither S t nor P e percolates and then 
both have divergent mean size (we shall call this third possibility in short 
formation of rings). 

Let us now briefly review our argument || that relates percolation prop- 
erties to the absence of a mass gap. Our statement was that if there was an 
equatorial strip S e that did not percolate for a certain c > 1 — e 2 /2, there could 
be no mass gap in the system. 

The argument is based on the imbedded Ising variables o~i = sgn(s(i)). 
Using these variables, the s.a. Hamiltonian becomes: 

= — cr^cr^jsii (j?) | - s ± (i) ■ s ± (j) (3) 

where sii(i) = s(i) ■ n and s±(i) = n x (s(i) x n). The c.a. model can be 
similarly described in terms of the variables cij, |sii(i)| and s±(i). In both 
models one thus obtains an induced Ising model for which the Fortuin-Kastleyn 
(FK) representation |L2| is applicable. In this representation the Ising system 
is mapped into a bond percolation problem: In the s.a. model a bond is 
placed between any like neighboring Ising spins with probability p = 1 — 
exp(— 2fls\\(i)s\\(j)). For the c.a. model a bond is also placed if after flipping 
one of the two neighboring Ising spins the constraint s(i) ■ s(j) > c is violated. 
From the FK representation is follows that the mean cluster size of the site 
clusters joined by occupied bonds (FK-clusters) is equal to the Ising magnetic 
susceptibility. In a massive phase the latter must remain finite. Hence, if 
the FK-clusters have divergent mean size, the original 0(3) ferromagnet must 
be massless (the Ising variables a are local functions of the originally spin 
variables s). 

Now notice that by construction for the c.a. model the FK-clusters with, 
say, a = +1 must contain all sites with s(i) ■ n > J(l — c)/2. Therefore the 
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model must be massless if clusters obeying this condition have divergent mean 
size. But the polar set P e consists of the two disjoint components P e + and P~ . 
For c > 1 — e 2 /2 there are no clusters containing elements of both P+ and P~ . 
Hence if for such values of c clusters of P e form rings, so do clusters of P e + 
separately and the 0(3) model must be massless by the argument just given. 

If we want to study the percolation or absence of percolation of the set A 
corresponding to a subset A C S N ^ of positive measure numerically, we have 
to consider a sequence of tori of increasing size L. On these tori we measure 
the mean cluster size of of A. If A percolates in the thermodynamic limit, by 
translation invariance (A) = 0(L 2 ); if its complement percolates (A) should 
approach a finite nonzero value, and if A forms rings we expect (A) = OIL 2 " 71 ) 
for some rj > 0. Therefore, if we define the ratio 

r=(P.)/(S e ), (4) 

for L — ► oo it should either go to if S e percolates or to oo if P e percolates; 
if neither S e nor P e percolates, then both form rings and the ratio r could 
diverge, go to or approach some finite, nonzero value depending upon the 
value of the critical index 77 for the two types of clusters. 

In the next section we will describe what our numerical simulations tell us 
about the percolation properties of equatorial strips and polar caps. 



3 Main numerical results 

Our results were obtained from a Monte Carlo (MC) investigation using an 
0(3) version of the Swendsen-Wang cluster algorithm |13| and consist of a 
minimum of 20,000 lattice configurations used for taking measurements. For 
each value of e we studied L = 160, 320 and 640 (for e = .78 we also studied 
L = 1280). 

In Fig.|l] we show the numerical value of the ratio r as function of c for 
/3 = for four values of e for the c.a. model on a T lattice. Three distinct 
regimes are manifest for each of the four values of e investigated: 

• For small c, r is increasing with L, presumably diverging to oo (region 
1). 

• For intermediate c, r is decreasing with L, presumably converging to 
(region 2). 

• For c sufficiently large depending upon e, r shows a very mild dependence 
upon L. 

This can only mean that for these values of e for small c P e percolates, 
for intermediate c S e percolates and for sufficiently large c both P e and S t 
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Figure 1: Ratio (P e )/(S e ) for various e values versus c 

form rings with quite similar (possibly equal) values of rj. Our data allow to 
deduce a semiquantitative 'phase diagram' in the (c, e)-plane of the percolation 
problem induced by the c.a. model on the T lattice for (3 = 0. 

This is shown in Fig.^. The solid line C is the curve c = 1 — e 2 /2; above 
that line the two polar caps cannot touch and therefore their union cannot 
percolate. The dashed line D represents the minimal equatorial width above 
which S e percolates. The point T at the intersection of the curves D and C 
gives an upper bound for c cr t, the value of c above which the c.a. model is 
massless. 

Let us explain how this picture was obtained: since for c = — 1 (no con- 
straint) the model reduces to independent site percolation, for which the per- 
colation threshold is known rigorously to be e = 1, curve D has to start at 
e = 1. With increasing c that threshold shifts to smaller values of e. Four 
points of the dashed line are in fact determined by the data displayed in Fig.|l|: 
the clearly identifiable four points where the lines for different lattice sizes 
L cross and the ratio r becomes independent of L determine the percolation 
threshold for the chosen value of e and so determine a point of the dashed line 
D. 

Two features of this diagram are worth emphasizing: 
1. An equatorial strip of width less than approximately e = .76 never perco- 
lates. 
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Figure 2: Phase diagram of the 0(3) model on the T lattice 

2. In Fig.[3] for approximately c > 0.4 a range of e's appears such that the 
ratio r again becomes approximately independent of L. This is signalling the 
appearance of a new 'phase' in which both S e and P e form rings (the dotted 
line separates it from the region of percolation of P e ). 

This regime of ring formation of both S e and P e is lying between the dotted 
and the dashed lines in Fig.||. Our data give strong evidence of its existence, 
but they do not determine in detail where the boundaries are. The dotted line 
has to run to e = for c = 1 because below it there is percolation of P e , and 
this is not possible above the solid line C e = e c (because it would conflict 
with the principle of non-coexistence of disjoint percolating sets). We drew 
the dashed line into upper right corner because we expect that eventually, for 
c approaching 1, any polar cap will start forming rings, thereby preventing 
percolation of the corresponding equatorial strips. 

But what it is essential for our conclusion that there is a massless phase is 
only that there is a regime below the dashed line D and above the solid line C, 
in which S e does not percolate and the two polar caps do not touch. In other 
words, the lines C and D have to cross (the crossing point is denoted by T in 
FigfJ). Since we found that for e < e = 0.76 S e never percolates, this means 
that for c > c eo = 0.71 the model is massless. In fact the massless phase must 
start earlier, and for instance based on our data we estimate that at c=0.61 
the c.a. model is already massless. 



7 



In the next section we will further corroborate the fact that for e < 0.76 
the equatorial strip does not percolate for any c. 

We would like to comment briefly on another recent paper dealing with 



percolation properties of equatorial strips in the 0(3) model: Alles et al f!4 



published a study showing that for e = 1.05 and (3 = 2.0 in the s.a. model S e 
percolates. Although strictly speaking our percolation argument applies only 
to the c.a. model, the result of Alles et al is not surprising since at (3 — 2.0 



the s.a. model is clearly in its massive phase [15], hence, by analogy with what 



happens in the c.a. model, one would expect that clusters of a sufficiently 
wide equatorial strip percolate (see [[RJ). The real issue, which the authors 
of JTJJ did not seem to appreciate, is whether in the c.a. model clusters of 
the equatorial strip S €c continue to percolate for c sufficiently close to 1. The 
numerics presented in Fig.|TJ suggest that that is not the case. 
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4 Corroborating numerics 



To corroborate our most important result, namely that for approximately e < 
.76 S e does not percolate for any value of c, we also measured (at j3 — 0) 
the ratio of the mean cluster size of the set P e t with e' = .5 to that of the 
set S e with e = 0.75 (e' was chosen so that P e t has equal density with S t ). 
The results are shown in (Fig.^J) . This figure shows that for c less than about 
0.4 (and greater than 0) the ratio grows very rapidly with L, indicating that 
P e t forms rings while S e has finite mean size; this region terminates around 
c = 0.4, where presumably also S e starts forming rings, and the dependence of 
the ratio upon L becomes much milder. Since for c > 0.4 the ratio continues 
to grow with L, at equal density, clusters of the polar cap are larger than 
those of the equatorial strip. The larger average cluster size of the polar cap 
compared to the strip of the same area is probably due to the fact that the 
strip has a larger boundary than the polar cap. This is in agreement with a 
general conjecture stated in ||, namely that for c sufficiently large, if two sets 
have equal area but different perimeters, the one with the smaller perimeter 
will eventually, for c approaching 1, have larger average cluster size. For the 
case at hand, this is apparently true for all values of c. 




Figure 3: The ratio of the mean cluster size of a polar cap of height .75 to that 
of an equatorial strip of the same height 
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5 Comparison to the 0(2) model 



The general belief, which we criticized in ref . J| , is that there is a fundamental 
difference between abelian and nonabelian models. To test this belief we also 
studied the ratio r for the c.a. 0(2) model on the T lattice. The phase diagram 
is shown in Fig.f|. Since in the 0(2) model the set V e can also be regarded 
as a set where e = a/4 — e 2 , certain features of that diagram follow from 
rigorous arguments. For instance it is clear that in the c.a. model there exist 
two intersecting curves C and C and in the region to their right the model 
must be massless ||, [7|. The precise location of the curves D (or D) must be 
determined numerically, something which we did not do. We did verify though 
that the ring formation region begins around c = —0.5. 




c 



Figure 4: The phase diagram for the 0(2) model on the T lattice 

6 Universality between the s.a. and c.a. mod- 
els 

In our opinion the arguments and numerical evidence provided so far give 
strong indications that the c.a. 0(3) model on the T lattice has a massless 
phase. Universality would suggest that a similar situation must exist for the 
s.a. models on the T and S lattices. To test universality we measured on 
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the S lattice the renormalized coupling both on thermodynamic lattices in the 
massive phase and in finite volume in the presumed critical regime (as in fl~7[l). 
Our data for the c.a. model on the S lattice only determine an interval (about 
.5 to .7) in which the massless phase of the model sets in; we tried to see 
if we could get a similar L dependence for the renormalized coupling in the 
s.a. model at a suitable (3 as for c = .61 in the c.a. model at (3 = 0. This 
seems to be indeed the case for (3 roughly 3.4. We went only up to L = 640, 
hence this equivalence between c and (3 should be considered only as a rough 
approximation, but there seems to be no doubt that the two models have the 
same continuum limit. 
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Figure 5: The step scaling curves of the s.a. and c.a. 0(3) models on the T 
lattice 



We also compared the step scaling curve in the s.a. and c.a. models. The 
step scaling curve is obtained as follows: on a periodic lattice of size L x L we 
define an apparent correlation length Namely let P = (p, 0), p = ^zl, 

n = 0, 1, 2, L — 1. Then define 

where 

G(p) = j- 2 (\s(P)\ 2 ); s(P)=J2e iP *s(x) (6) 
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Leaving f3 respectively c fixed, one doubles L and measures also £(2L). The 
step scaling function gives the ratio 2£(L)/£(2L) versus L/£(L). In the con- 
tinuum limit [L — > oo and /3 — > /3 cri at L/£(L) fixed), this procedure produces 
a unique curve characterizing the universality class of the model. In Fig.|5| 
we present the step scaling function for the c.a. and s.a. models. The data 
were produced by adjusting (3 and c so that at L = 20 we obtain roughly 
the same £(L) in the two models. After that, leaving (3 respectively c fixed, 
L was doubled until L = 320. As can be seen, the two step scaling curves 
agree reasonably well; the slight disagreement is probably due to the fact that 
the two curves have to agree only in the continuum limit (£ — > oo), i.e. they 
have different lattice artefacts, whereas the largest value of £ reached was only 
approximately 35 lattice units. 



7 Heuristic explanation of the transition 

It is intersting to note that there is a heuristic explanation for both the ex- 
istence of a massless phase in the s.a. 0(3) model and for the value of /3 crt . 
Indeed it is known rigorously that in 2D a continuous symmetry cannot be 
broken at any finite f3. In a previuos paper || we showed that the dominant 
configurations at large (3 are not instantons but superinstantons (s.i.). In prin- 
ciple both instantons and s.i. could enforce the 0(3) symmetry. In a box of 
diameter R the former have a minimal energy E inst = 4ir [Tj|] while the latter 
E s .i. — 5 2 7r/ lni?, where 5 is the angle by which the spin has rotated over the 
distance R. Now suppose that (3 crt is sufficiently large for classical configura- 
tions to be dominant. Then let us choose 5 = 2ir (restoration of symmetry) 
and ask how large must R be so that the superinstanton configuration has 
the same energy as one instanton. One finds \nR = tt 2 . But in the Gaussian 
approximation 

(s(O)-s(x)) Inx (7) 

jJTT 

Thus restoration of symmetry occurs for lnx ~ ir/3. This simpleminded 
argument suggests that for (3 > tt instantons become less important than s.i.. 
Now in a gas of s.i. the image of any small patch of the sphere forms rings, 
hence the system is massless. While this is not a quantitative argument, we 
believe it captures qualitatively what happens: a transition from localized 
defects (instantons) to super-instantons. 
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8 Absence of asymptotic freedom 



Next let us discuss the connection between a finite (3 crt and the absence of 
asymptotic freedom. It follows from our earlier work concerning the conformal 
properties of the critical 0(2) model [1C]. We refer the reader for details to that 
paper and give only an outline of the argument. The s.a. lattice 0(N) model 
possesses a conserved isospin current. This currrent can be decomposed into a 
transverse and longitudinal part. Let F T (p) and F L (p) denote the thermody- 
namic values of the 2-point functions of the transverse and longitudinal parts 
at momentum p, respectively. Using reflection positivity and a Ward identity 
we proved that in the massive continuum limit the following inequalities must 
hold for any p ^ 0: 

< F T (p) < F T (0) = F L (0) < F L (p) = 2(3E/N (8) 

Here E is the expectation value of the energy 

E = (s(i) • s(j)) 

at inverse temperature (3. Since E < 1 it follws that if f3 crt < oo F T (0) — F T (p) 
cannot diverge for p — > oo as required by perturbative asymptotic freedom |1| . 
In fact, for (3 crt = 3.4 (which is a reasonable guess) F T (p) must be less than 
2.27, in violation of the form factor computation giving F T (0)—F T (oo) > 3.651 

0- 



9 Concluding remarks 

Since the implications of our result, that for the c.a. model a sufficiently narrow 
equatorial strip never percolates, are so dramatic, the reader may wonder how 
credible are the numerics. The only debatable point is whether our results 
represent the true thermodynamic behaviour for L — > oo or are merely small 
volume artefacts. While we cannot rule out rigorously the latter possibility, 
certain features of the data make it highly implausible: 

• Small volume effects should set in gradually, while the data in Fig.|l| 
indicate a rather abrupt change from a region where r is decreasing with 
L to one where r is essentially independent of L. 

• For c — > 1 at fixed L, r must approach the 'geometric' value r — 2/e — l. 
As can be seen, in all the cases studied, throughout the 'ring' region r is 
clearly larger than this value, while it should go to if S e percolated. 

• In Fig.|3] there is no indication of the ratio going to for increasing L. 
Moreover the dramatic change in slope around c = .4 indicates that 
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the polar cap P e > starts forming rings at a smaller value of c than the 
equatorial strip S e . 




3 4 5 6 7 

ln(L) 



Figure 6: Mean size of clusters of the polar cap P + .1 at c = — 1 (Bernoulli) 
and at c = versus L. The correlation length £ is approximately 53 lattice 
units. 

We have additional numerical evidence that clusters of a polar cap P e + 
smaller than a hemisphere (s ■ n > e/2 > 0) form rings for some c < 1. 
Namely we investigated the case e = 0.1. For the the case c = — 1 (Bernoulli 
percolation) it is known rigorously that clusters of this set have finite mean 
size. As can be seen from Fig|| our numerical values at c = — 1 corroborate this 
fact. In the same figure we show the mean cluster size of cluster of P + x at c = 0, 
where the correlation length is approximately 53 lattice units. Even though 
we increased L up to 1280, the mean cluster size shows no sign of leveling off, 
growing in fact like some power of L, consistent with the formation of rings. 

Therefore there is good numerical evidence that for c = 0, where we can 
reach the thermodynamic limit, clusters of this polar cap form rings. The 
natural expectation would be that the mean cluster size of a subset of a hemi- 
sphere is a nondecreasing function of c. This is borne out by the numerics, 
as shown in Fig.|7|. There we represent the mean cluster size of P + i a ^ fixed 
L = 640 function of c. The data support the assertion that for any c > 
the mean size of the clusters of P + x diverges, which, via our argument, implies 
that the c.a. model must cease being massive for some c < 1. 

Thus we doubt very much that the effects we are seeing represent small 
volume artefacts. Moreover, if s z , the ^-component of the spin s remained 
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Figure 7: Mean size of clusters of the polar cap Pq~.1 at L — 640 versus c. 

massive at low temperature and in fact an arbitrarily narrow equatorial strip 
percolated, one would have to explain away our old paradox @, 0: if such a 
narrow strip percolated, an even larger strip would percolate and on it one 
would have an induced 0(2) model in its massless phase, in contradiction to 
the Mermin- Wagner theorem. 

Consequently it seems unavoidable to conclude that the phase diagram 
in Fig.|2| represents the truth, that a soft phase exists both in the s.a. and 
the c.a. model and that the massive continuum limit of the 0(3) model is 
not asymptotically free. In a previous paper || we have already shown that 
in nonabelian models even at short distances perturbation theory produces 
ambiguous answers. The present result sharpens that result by eliminating 
the possibility of asymptotic freedom in the massive continuum limit. 

Acknowledgement: AP is grateful to the Humboldt Foundation for a Senior 
US Scientist Award and to the Werner-Heisenberg-Institut for its hospitality. 
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